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ABSTRACT 



This investigation analyzed families of ascent 
trajectories to determine total velocity .require- 
ments, division of velocity impulses and the ex- 
pected errors in position and velocity upon arrival 
in the terminal region caused by uncertainties in 
the position, magnitude and direction of the 
velocity vector at burnout of the booster stages. 

Conditions at burnout of the booster and apogee 
of the trajectory were solved using an IBM 650 
computer; the error analysis was accomplished by 
numerical integration on IBM 709 equipment repeated 
up to the apogee time of the zero-error trajectory. 

From the tabulated and graphical results pre- 
sented herein, some trade-offs between propulsion 
requirements and the accuracy of guidance and con- 
trol systems can be seen. 
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OBJECT 

I 



It is the object of this thesis to 
present in graphical form the effects 
of errors in position and velocity at 
booster burnout on the position and 
velocity of a rendezvous vehicle at 

i 

zero-error apogee time. 
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CHAPTER I 
INTRODUCTION 



Requirements 

Long before the launching of the first earth 
satellite the military, scientific, economic and 
political ramifications of such a vehicle were con- 
sidered and discussed. Each ©f the above considera- 
tions generates adequate reasons for a satellite 
rendezvous . 

The most Important and probably the most 
immediate requirement is in the field of the mili- 
tary. In orbit about the earth today are satellites 
capable of reconnaissance; little additional equip- 
ment is required to extend this capability to ferret 
and jamming missions or ICBM early warning. Inter- 
rogation and possible negation of such enemy vehicles 
are missions that require a satellite rendezvous. 

Satellite stations in space require the physical 
exchange of personnel, equipment and fuel. Com- 
pletion of such an exchange requires a "soft" rendezvous 
and continuation of the two vehicles in the same orbit 
until the transfer is complete. 
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Both civil and military organizations are 
developing communications satellites; these de- 
vices involve much complex and expensive equip- 
ment that will, in time, require repair or re- 
placement. Satellite rendezvous is a solution 
to this problem. 

Communications satellites are given special 
consideration in this report by the inclusion of 
data for orbit altitudes as great as 20,000 nauti- 
cal miles above the surface of the earth. Earth 
satellites much above this altitude are of little 
interest since they are subject to large pertur- 
bations by the sun and moon and hold no advantage 
over the satellite with a period of 2 4 hours. The 
effect of such perturbations and the difficulties 
to be encountered in maintaining a high altitude 
orbit are discussed by Edelbaum in reference 1. 

Background 

The subject of satellite rendezvous has re- 
ceived its share of attention in the past several 
years. Literally hundreds of papers and articles 
have been written on the various phases of rendez- 
vous and interception. Most of these works are 
very general and deal primarily with ideas rather 
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than figures and hardware. Other authors look into 
the problem in detail and find it to be enormous; 
it seems to take a well-funded project to generate 
any amount of quantitative date on a particular 
phase. Many studies are initiated on the basis 
that the proposed payload and vehicle character- 
istics are known and the task is to find out how 
(or if) the desired mission can be accomplished. 
Much of this work is security classified and is, 
therefore, unaccessible. It is the purpose of 
this thesis to conduct an academic yet practical 
study of satellite rendezvous so that students may 
gain an insight into a very complex problem. 

An important part of this subject that is 
specifically not included in this thesis is the 
topic of guidance and control during the terminal 
phase of the rendezvous maneuver. This topic 
appears to be a study in itself inasmuch as the 
terminal phase shifts to "closed loop" control 
with feedback from radar as the rendezvous vehicle 
closes on the target. Some facts about the capa- 
bilities of the acquisition system and, in fact, 
entire specification of the means to be used by 
the rendezvous vehicle in completing the terminal 
phase must be outlined to justify the definition 
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of a successful arrival. Sears and Felleman in re- 
ference 2 give an excellent example of this type 
study. Nason establishes the conceptual limitation 
of proportional navigation in reference 3. 
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CHAPTER II 



STATEMENT OP THE PROBLEM 



General 

The earth satellite Is a unique device 
since , unlike vehicles that operate In the earth's 
atmosphere, its flight is a regular pattern in 
space determined primarily by the velocity imparted 
to the vehicle at the beginning of its orbit. While 
velocities of earth satellites are high compared to 
any motion within the earth's atmosphere, the motion 
of the satellite is predicted with accuracy and it 
Is this statement that forms an important base uti- 
lized In this thesis. 

Satellite rendezvous is the problem of bringing 
together two objects of relatively small size when 
the original distance between the two objects and 
their velocities are both very large by ordinary 
engineering standards. When the target orbit ele- 
ments are known, the task of the rendezvous vehicle 
(RV) Is to match them. In practically any clrcum- 
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stances, even the simplest case, there exist an 
infinite number of possibilities for accomplish- 
ing a rendezvous . 

The general case of matching the ephemeris 
of an arbitrary target after launch from an 
arbitrary point in space (surface of the earth) 
is unmanageable and probably unrealistic. Al- 
though an effort is made to keep the problem as 
general as possible, certain conditions are set 
before proceeding. 

The Specific Problem 

The study begins with the assumption that the 
orbit of the target and the trajectory of the RV 
are in the same plane. If the orbit of the target 
has an Inclination less than the latitude of the 
launching point then plane -change Is required as 
the target never passes over the launch point. 

With angles of inclination greater than this lati- 
tude the launch point is in the plane of the target 
twice each day giving two opportunities for launch 
into the target plane. Except for polar orbits, 
one of these opportunities requires a westward 
launch, an undesirable choice, thus the desirable 
eastward launch into the target plane is available 
only once a day. 



Now Consider further the odds that the target 
is at a suitable position for rendezvous at the . 
time the launch point passes through the target 
orbit. Obviously a wait of many days may be re- 
quired before the correct conditions are met . 

A study of the wait time problem is made by 
Houbolt in reference 4. In this paper he de- 
monstrates for a specific case the limited number 
of rendezvous opportunities available in a given 
period and points out the increase in this number 
given by the "Adjacent Launching Technique". In 
this scheme the RV is launched into an orbit of 
the same altitude as the target with the target 
at a suitable position for rendezvous. When the 
RV crosses the target orbit it adapts to this 
new orbit with a plane-change maneuver and begins 
the terminal phase of the rendezvous. 

Inasmuch as changes of orbit inclination are 
accomplished by an impulse normal to the original 
velocity vector, the change of plane requirement 
Is in addition to the in-plane requirements and 
is not considered further. 

Since the target is, by definition, the point 
in space to which we wish to move the RV, it Is 
reasonable- to assume that the target has a regular 
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motion easily predicted. A circular orbit best 
meets this requirement and this is not an un- 
realistic assumption; for many reasons an approxi- 
mately circular orbit is desired for reconnaissance, 
survey or communications satellites. 

Now within a given plane and with a circular 
target orbit there still remain an unlimited 
number of trajectories that meet the requirement 
of placing the RV at the satellite's position in 
3pace. These trajectories form the basic conic 
sections, ellipse, parabola and hyperbola, passing 
through the two required points. Only a minimum 
of computation is required to see that hyperbolic 
and parabolic trajectories have excessive velocity 
requirements when compared with the elliptical 
path. The elliptical trajectories that qualify are 
those with apogee greater than or equal to the al- 
titude of the target. The former, often called a 
"fast ellipse", crosses the target orbit and this 
is undesirable for two reasons? 

1. To match the orbit of the target the KV 
requires in -plane normal thrusting at the cross- 
over point. Normal thrusting is inefficient and 
is only used as a last resort. 
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2. Closing velocities at the cross-over are 
very high and rendezvous requires critical timing. 

The ascent ellipses considered in this thesis 
are those with apogee equal to the altitude of the 
target. This type trajectory gives a flight path 
at apogee that is tangent to the orbit of the 
target . 

Within the conditions of being in-plane and 
tangent at arrival/ there is an alternative, four- 
impulse, "parking orbit" ascent. Briefly, parking 
orbit ascents consist of a launch into a low alti- 
tude orbit, a waiting period until the relative 
positions of RV and target are correct, then a 
transfer to the target orbit. The primary advantage 
of such a scheme is the relaxation of strict re- 
quirements for leaving the launching pad at a given 
time; the low altitude orbit allows the RV to "gain" 
on the target. The disadvantages of such a system 
are related to this thesis in that errors in burnout 
velocity are not only multiplied by the requirement 
for two "transfer" maneuvers, but initial errors in 
many cases "grow" during the wait period. These 
and other factors are considered in detail by Duke, 
Goldberg and Pfeffer in reference 5 . 
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Coasting In free flight from burnout of the 
booster stages until the KV is at the altitude 
of the target is termed a "direct ascent". Such 
a trajectory has a minimum velocity requirement 
in the classical Hohmann maneuver and an upper 
limit in the vertical ascent. This thesis will 
be concerned with these two-impulse type rendez- 
vous trajectories. 

To review, then, the specific problem in- 
vestigated in this thesis is an in-plane rendezvous 
with a target in circular orbit utilizing a direct 
ascent elliptical trajectory that is tangent to the 
target orbit. 

Basic Assumptions 

The earth is assumed to be a homogeneous 
sphere. The earth is, of course, quite irregular; 
it is more an ellipsoid than a sphere. 

This oblateness has three main perturbatlng 
effects upon the orbit of an earth satellite, all 
dependent upon altitude and inclination of the 
orbit : 

1. Regression of the nodes. This is motion 
of the orbital plane about the earth's axis in a 
direction opposite to that of the satellite. 
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2. Precession of the line of apsides. This 
motion is an angular rotation of the axis of the 
orbit in inertial space. It may be in the same 
or the opposite direction of the satellite motion 
depending upon the amount of Inclination. 

3. The oblateness of the earth also causes 
slight fluctuations of the radius (altitude) and 
period of the satellite. 

All of the above effects are important in 
sustained satellite operations and their effects 
on orbits adds to the knowledge of the shape of 
the earth, however, such effects are of small 
magnitude in direct ascent rendezvous trajectories. 

As an example, in reference 6 the formula: 

, R 2 

A - 2.t ( 1.637 x 10"^) cos i 

( R e + Z s ) 2 

is given as the measure of angular motion of the 
nodes in radians per revolution of the satellite. 
For an altitude of 1000 n. miles and with a cir- 
cular orbit inclined at an arbitrary angle of 45°, 
the regression of the nodes is 4.38 x 10“^ radians 
for every revolution of the target. The ascent 
trajectories described in this thesis are of short 
duration compared with the period of the satellite 
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target and only a portion of this regression of 
the nodes takes place while the KV is In flight. 

For the 1000 n. mile target the ascent time varies 
from 3154 seconds to 843 seconds or 11# to 40# of 
the period of the target. Thus the regression of 
the nodes is from 13,000 to 47,000 feet at this 
geocentric radius. 

This example Is given only to show the order 
of magnitude of the motion in this case. In the 
rendezvous situation the perturbations affect 
the HV as well as the target at the same time and 
in the same direction. The error is, then, of 
differential order when compared with the above. 

1 

If these effects are expressed as functions 
of position they can be introduced into the inte- 
gration of the equations of motion for the 
rendezvous computations, but they are not the sub- 
ject of the present investigation and can influence 
it only through higher order effects. 

To continue with the basic assumptions, 
booster burnout altitude is taken as sixty nautical 
miles for all cases. Without a detailed examination 
of the booster performance and maneuvering scheme 
it cannot be determined whether this altitude is 
appropriate for a given trajectory, but some arbi- 
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trary altitude is needed for a common base from 
which to compare the coasting performance. This 
altitude is assumed in numerous references such 
as Eggleston in reference 7- 

It is further assumed that the ephemeris of 
the target is known exactly. The technique of 
tracking an earth satellite and determining its 
ephemeris is fairly well formalized and high 

-i l 

speed digital computer routines are available for 
this purpose. The technique ifl described by Saindorff 
in reference 8. For this study, however, uncertainties 
of the satellite position are simply uncertainties in 
the origin of the final error coordinates. While this 
may be an important consideration, it is independent 
and can be introduced as necessary in a system study. 

Although the ability to predict launch time 
continually improves . with experience in missile 
operations, there are still uncertainties in time 
of lift-off as well as error build-up during burn- 
ing. The time of flight or coasting time varies 
considerably with the particular trajectory 
selected (as seen below) giving some flexibility 
in matching time. A pre-launch computer keeps 
track of this time and, over some useful range 



determined by other considerations, selects the 
exact trajectory to match time. In a similar 
manner, during the burning period performance is 
monitored and cut-off corrections computed to 
match the time. In very general terms, Duke, 
Goldberg and Pfeffer in reference 5 take notice 
of this type error but its implementation de- 
pends upon system specifications. If an error 
in time alone is accepted and no compensation 
is made then the arrival error is easily com- 
puted from the satellite velocity. To get to 
the point of the study, it is assumed further 
that booster burnout is compelled to occur ex- 
actly on time. 

Uncertainties Considered 

The uncertainties considered in this thesis 
are those in the booster burnout velocity vector 
These errors are described as; 

1. Position errors 

a. ahead (RANGE error) 

b. above (ALTITUDE error) 

c. to one side (TRACK error) 

2. Velocity errors 



a. etfror in magnitude 
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b. error in direction 

(1) deflected above or below 
the ideal 

(2) deflected to one side of 
the ideal 

The exact meaning of "position errors" de- 
pends upon the type of guidance equipment utilized 
by the booster stage; the full implication of this 
point is carefully discussed in Chapter III. 
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CHAPTER III 

COMPUTATIONAL METHODS AND PROCEDURES 

Analytical Computation of Reference Trajectories 
The elements that describe the target orbit 
and the family of possible direct ascent ellipses 
tangent to it are defined with the help of Fig. 1 
below. 




Fig. 1. Definition of Elements Used in Computational 
Procedures. 
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Prom the figure it is seen that selecting 
a target altitude fixes R<* for all ascent ellipses, 
since the rendezvous condition is that apogee 
occurs at target altitude. An additional variable 
is needed to specify a particular trajectory from 
the possible family of ellipses, and the range 
angle, <$> , best describes what is a long or short 
ascent. Although these two variables are selected 
for classification of results, the conventional 
two-body equations are more familiar in terms of 
the standard parameters e, a and 0. The desired 
values are computed by successive substitution in- 
to the numbered equations below. 

These equations are not involved, but they 
require careful numerical evaluation for accurate 
results. With dimensions in millions of feet, 
the IBM 650 carrying eight significant digits is 
adequate for computation directly in engineering 
units. The three physical constants for the 
numerical work are: 

G (gravitational constant times mass 
111 of the earth) 

1.40768 x 10 l6 ft 5 /sec 2 Ref. 8 

R (radius of a sphere with volume 
e equivalent to that of the earth) 

3437.866 nautical miles 



Ref. 9 



International nautical mile * j 6076.1155 feet 

Ref. 10 
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Altitude traversed in ascent 

z * - z s - z 1 
Minimum eccentricity 

e « z * 
min 2 R- l + z* 

Angle from perigee , 

9 - (180° - 4> ) 



(3.1) 



(3.2) 



(3.3) 



Eccentricity 



e » — 



cosO +1 \ 


1_ 


2 

' “N 

cosQ + 1 


4- 


z */ r i "I 


2 [ cosS+1+ \ 


4 


COS9+1+ t=> — 

r i 


t 


„ n z * 
cos9+l+ jp~ 

^ 1 > 



Satellite velocity 



Gm 



R 1 +z * 



Booster burnout velocity 




Apogee velocity 



+ 1-e 



(3.4) 



(3-5) 



(3.6) 



(3.7) 
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Injection velocity increment 



V - V, 
s 



Semi -major axis 



% + Z oC 



" 1 + e 

Booster burnout attitude angle 



' 1 $ - arctan 
Time to apogee 



e sinQ 
1 + e cos© 



■Nfct' - arccos ( — — - ) + e sin ^ arccos ( 



(3.8) 



(3.9) 



( 3 . 10 ) 



a - R, 



ae 
(3.11) 



This information, computed for 12 target alti- 
tudes from 200 to 20,000 nautical miles and range 
angles of 180, 90, 45 and 22.5 degrees permits 
selection of a smaller representative sample for 
detailed error analysis. Target altitudes of 400, 
1,000, 5 ,000 and 20,000 n. miles are used in the 
error analysis to follow. 

The Error Analysis 

A complete error analysis of even this basic 
set of trajectories by analytical methods, if 
possible at all, is certainly complex because 
several orbital elements are often changed by a 
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single error at booster burnout. For instance, 
the variation in semi -major axis can be found as 
a function of a change in If , but before an error 
can be isolated, the corresponding change in ec- 
centricity and the shift of the line of apsides 
must be determined. Finally, arbitrary position 
as a function of time is trancendental. 

Braham and Skidmore in Ref. 11 give analytical 
expressions for changes in altitude, final velocity 
and final flight path angle as a function of certain 

initial condition errors. These expressions are 

* < 

cumbersome. The one for final flight path angle 
requires evaluating six lengthy expressions to get 
the necessary partial derivatives, which in turn 
require evaluating four other equations equally 
as long to solve for some of the coefficients. 

Even if a complete set were available, evaluation 
in any quantity calls for a machine computer pro- 
gram. 

Such analytical solutions are interesting 
from the theoretical point of view, but it is 
felt that the most practical way to get a com- 
plete error analysis is by numerical integration 
of the equations of motion. 
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The motion of a body subject to an accelera- 
tion is described by the vector equations: 



V - 
R - 



Ja at + V 0 



J V dt + R 



o 



( 3 . 12 ) 

(3.13) 



In coasting flight near the earth, the only 
acceleration is that due to the gravitational 
field, G, which is a function of R. Thus, for 
step-by-step integration, the two simultaneous 
equations are actually solved in order as follows: 



a - 




(3.14) 



This acceleration is expressed in components in 
whatever coordinates V is given, and each com- 
ponent of V and R is computed as: 




R 



i 




At + V ±0 


(3.15) 


At + R l0 


( 3 . 16 ) 



The resulting R is used to compute a new G; the 
resulting V is the V Q for the next step, and the 
resulting R is also the R Q for the next step. This 
procedure is repeated up to apogee time. 

A complete program with this capability, de- 
veloped by Mr. Robert Brown of the M.I.T. InBtru- 
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mentation Laboratory and adapted to this problem, 
is used in the error analysis. With this program 
and the IBM 709 computer at the M.I.T. Computation 
Center, trajectories for four altitudes and four 
range angles are thoroughly investigated. 

The method is straight forward and accurate 
when used with a suitable integrating time step. 

The validity can be checked for each trajectory 

- 

by integrating over the complete time of flight 
with no errors introduced and comparing the 
apogee conditions with the analytical results 
from the equations above. As an example of the 
accuracy attained the zero-error numerical in- 
tegration for ascent to a 5>000 n. mile altitude 
at a range angle of 90° using a 50 second time 
step agreed to within 80 feet ALTITUDE and 7 feet 
RANGE distance - this is at a geocentric radius 
of over 51 million feet. 

The program uses spherical coordinates in 
the actual integration, but keeps track of the 
components in terms of an orthogonal frame called 
a platform frame. This frame is introduced to 
simplify the presentation of the error components, 
and has a positive X axis in the RANGE direction 
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along the horizontal in the direction of flight, 
a positive Z axis upward along the vertical 
ALTITUDE direction, and a positive Y axis to the 
right looking in the direction of flight. 

A platform frame is specified because for 
inertially guided vehicles directions are re- 
ferred to a stabilized platform. The platform 
is non-rotating in earth-centered inertial co- 
ordinates (type l) or its Z axis tracks the local 
vertical (type 2). The effect of position error 
on both types is shown in Pig. 2. 

The first type seems most likely, for this 
application, but it is used here for a more basic 
reason. With error effects available for the 
first type, they can be computed for the second 
type. Note from the figure that ALTITUDE errors 
alone are the same in either system, but that 
with RANGE or TRACK error, the relative alignment 
of the G vector depends on the type system assumed. 
It is seen though, that a type 2 RANGE error con- 
sists of a type 1 RANGE error plus a velocity mis- 
alignment ( ) which is easily computed from 

booster burnout radius and the RANGE error. 

Effects due to the actual change in height 
caused by a type 1 RANGE or TRACK error are negli- 
gible because, for reasonable errors, this change 



T 
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in height is very small, and the corresponding 
position error derivatives are also very small. 
The same is not true for misalignment errors. 

The errors may be fractions of a milliradian, 
but this kind of error has a large end effect. 

Finally, for TRACK error, the motion due 
to the type 1 position error must be obtained 
by numerical integration, but for the type 2 
platform where the ALTITUDE component of velocity 
is aligned with the local vertical, the entire 
motion lies in the plane of the error trajectory. 
When a trajectory lies entirely in a plane and 
the inclination with respect to another plane 
is known, the lateral displacement is just a 
function of the geometry. This point is brought 
out in more detail in Chapter IV. 
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Comparison of two types of platform coordinate 
frames. Angle for reasonable position error 
greatly exaggerated. 




Ref. local horizontal 




Error in magnitude of 
booster burnout velocity. 




Track Error - Type 1 
system as used herein. 



In-plane misalignment 
of booster burnout 
velocity . 




Lateral Misalignment 
of booster burnout 
velocities 



Pig. 2. Platform Coordinates and Definitions of Booster 
Burnout Errors . 
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CHAPTER IV 
RESULTS 



Summation of Velocities 

The summation of booster burnout velocity and 
Injection velocity Increment (XV) for various alti- 
tudes and range angles is shown in Fig. 3. It must 
be borne in mind that the values represented here 
are the theoretical minima; they do not reflect the 
additional velocity that must be provided to over- 
come drag and gravity losses during accelerations 
to V-^. 

What is considered important about these plots 
is the effect of altitude and range angle on XV. 
Note that while analysis always prescribes the 
Hohmann maneuver for minimum XV, this is not 
necessarily the primary consideration in a practical 
case. For instance, with a 400 n. mile orbit alti- 
tude the range angle is reduced to 55° before XV 
Increases 1# over the Hohmann value; it is reduced 
to 18 ° before XV increases 5$. Similarly for a 
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Summation of Booster Burnout Velocity and Injection Velocity 
Increment. Altitudes in nautical miles are shown adjacent 
to the curves. 



29 



1000 n. mile orbit altitude the range angle is 
down to 40° before XV increases by 5$. 

These examples are cited only to call at- 
tention to what may be one corner of a many 
sided compromise. 

Velocity Division 

Pig. 4 is a plot of the separate components 
of XV versus range angle for a number of orbit 
altitudes. The information is displayed on two 
separate plots for low and high altitudes to 
eliminate the clutter of several overlaps. 

V^, the booster burnout velocity, is the 
velocity which must be imparted to the RV to 
send it coasting up to apogee. V^, the injection 
velocity increment, represents that propulsive 
capability that must be carried to the rendezvous 
point in order that satellite velocity be matched. 

Prom another point of view, V^ is the relative 
velocity between target and RV. If the ascent 
trajectory is perfect and time matched exactly, 
the injection thrusting can be initiated by elapsed 
time, but with error buildup during the ascent it 
seems most likely that the terminal phase will be 





^* Effect of Altitude and Range Angle on BooBter 
Burnout Velocity and Injection Velocity 
Increment . 



51 



closed -loop. Thus, this relative velocity may 
be important during the acquisition olf the target 
by the HV. 

As before, for computation of propulsive re- 
quirements the drag and gravitation losses must 
be added to Then, once an effective exhaust 

velocity or specific impulse is known or assumed, 
quick preliminary calculations can be made con- 
cerning vehicle design. These might consist of 
selecting an altitude and/or range angle within 
the capability of a certain system, or, in deter- 
mining what configuration and staging or what re- 
start capability is needed to accomplish a given 
mission. 

This information forms the second part of 
the propulsive picture, and another side of the 
overall compromise involved. 

Error Derivatives 

The results of the error analysis are re- 
duced to partial derivatives since it is shown 
that the variation of a final component with re- 
spect to a change in the initial value of any 
other component is linear over what is considered 
to be a generous range of variation. This direct 
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proportionality is verified for the 1000 mile 
orbit for speed errors up to 60 feet per 
second, for angular deviations of the velocity 
up to two milliradians and for position errors 
up to 12,000 feet in any direction. As shown 
for the errors in booster burnout velocity, 
plotted on Pig. 5, there is no tendency to de- 
part from linearity at 60 feet per second. Be- 
cause the functions being considered are smooth 
and continuous, they do not bend abruptly so it 
is considered that booster burnout velocity 
errors up to 100 feet per second and angular 
deviations involving components of that size 
could be entertained with confidence. It does 
not seem likely that such large errors are 
realistic, although what is tolerable certainly 
depends on what is said about the terminal 
system capabilities. In any case, the additional 
computer time required to investigate by this 
method and establish the absolute limit of appli- 
cability is not justifiable with the many other 
items being studied with the computer. 

The types of booster burnout uncertainties 
described above are classifiable as in-plane and 
lateral errors. With three components each to 



Range Velocity Error, ft/sec X Range Error 
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Fig. 5. Effect of Errors in the Magnitude of the 
Booster Burnout Velocity on Range and 
Range Velocity Demonstrating Linearity. 
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define the position and velocity vectors, six 
differentials are sufficient to describe any 
possible error in the booster burnout condition. 

To express the reference apogee or final value 
of each of the six vector components in terms 
of each of the six initial compohents leads to 
56 partial derivatives. There is, however, 
clearly no way for the in-plane effects to in- 
duce any lateral perturbations and, as expected, 
the lateral errors introduce only infinitesimal 
effects on the final in-plane performance. With 
this independence, the number of error derivatives 
is reduced to 20. 

While it seems more natural to describe the 
Initial velocity errors in terms of variations 
in magnitude and direction, the final velocity 
errors are given in terms of orthogonal components. 
This is done primarily for ease of combination of 
components as demonstrated later. Also, since both 
the target and the reference apogee point have zero 
Y and Z components of velocity, the computed error 
components are the actual relative velocity in the 
ALTITUDE and TRACK directions as viewed from the 
satellite. 

The X component of the error velocity, while 
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generally quite small compared to either 
the actual or relative velocities in the RANGE 
direction, is still, important . It is a part of 
the terminal problem, and the rendezvous is not 
complete until all components of relative 
velocity have been exactly matched. 

The derivatives relating the initial and 
final variations are made up of combinations 
of the following: 

1. In-plane - Initial Position (x and z) 
and velocity (v^ and T ) 

Final position (X and Z) add velocity 
(V x and V z ) . 

2. Lateral - Initial Position (y) and velocity ( f ) . 
Final position (Y) and velocity (Vy). 

The general form of the drivative as used for 
evaluating ascent errors is y f ^ j - This 

does not preclude using the inverse.' If, for 
instance, it is desired to define allowable final 
error first, then th© permissible booster burnout 
deviations of a type could be found. However, 
there are, in general, many combinations of booster 
burnout errors that can combine to cause similar 



final errors. 



The derivatives are plotted in Figs. 6 through 
15 with range angle as the primary variable 
and the target altitude as the step parameter. 
Each is explicitly defined on the respective 
plot with units in feet, feet per second and 
milliradians . 

It is important to note that as range 
angle approaches zero, the eccentricity of the 
trajectory ellipse approaches unity and the 

regular forms of the equations become very 

/ 

sensitive. Judging from a few of the cases 
for which the vertical ascent may be evaluated, 
it appears that the error functions in fact are 
quite unstable at low range angles. For this 
reason, none of the results should be extra- 
polated below the 22.5° range angle where the 
error computations were cut off. There seem 
to be practical reasons for avoiding very steep 
ascents such as the high relative velocities, 
highest possible 2LV, and need for rapid change 
of attitude in the terminal phase. For such 
reasons, it is considered that the lack of 
coverage at these short angles is not important. 
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Pig. 6. Effects of Range Error at Booster Burnout on Position at Zero-error Apogee Time 
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Fig. 8. Effects of Altitude Error at Booster Burnout on Position at Zero-error Apogee Time 
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Fig. 10. Effects of Error in Magnitude of 3ooster Burnout Velocity on Position at 
Zero-error Apogee Time. 
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Fig. 15. Effects of Lateral Misalignment of Booster Burnout Velocity on Position 
Velocity at Zero-error Apogee Time. 
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Independence of Derivatives 

A final important result is the independence of 
the in-plane derivatives. To investigate all possible 
combinations of booster burnout errors in varying 
magnitude involves endless computation, but the total 
error is shotftt in Table I to be made up of the sum of 
the errors due to the separate effects. Table I is 
compiled by introducing a random combination of errors 
into the 1000 n. mile, ’ 90 ° range angle trajectory and 
checking the computed output for component addition 
from the separate derivatives. 



TABLE I 

Pinal Error Components 



Booster Burnout 
Error 


X 

feet 


Z 

feet 


V x 

ft/sec 


V z 

ft/sec 


(l) +20 ft/sec 


-44,293 


58,693 


-79.94 


40.25 


(2) -2 mr 


-103,095 


40,404 


-73.15 


-10.03 


(3) 1 mile ahead 


-7,711 


6,048 


-10.92 


-.75 


Sum of (1) + (2) 


-147,388 


99,097 


-153.09 


30.22 


Computed Results 


-147,387 


98,929 


-153.10 


29.89 


Sum of (1) + (3) 


-52,004 


64,741 


- 90.86 


39.50 


Computed Results 


-51,994 


64,753 


-90.87 


39.48 


Sum of 

(1) + (2) + (3) 


-155,099 


105,145 


-164.01 


29.47 


Computed Results 


-155,080 


104,964 


-164.04 


29.09 
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Let the assumed booster burnout errors be de- 
signated by £x, 8 z, h'V 1 and &Y for in-plane 
errors and by S y and S for lateral errors . Now 

if the final errors are designated by AX, AZ, AV 

A 

and AV Z for in -plane errors and by AY and AV.^ for 
lateral errors then a complete final error for a 
booster burnout error in, say, speed, range and 
lateral alignment is: 



AX 



ax 

3V-, 



* V 1 





S'Vi 


+ _az 

ax 


Sx 


BV„ 

AV - 

x dV 1 


Sv x 


3V x 

+ 

ax 


8x 



3 V Z c 

-377 & v i + 



av, 



si 



AY 



AV. 



dC 0 

iic* 

d<r 



where the partial derivatives are read from the 
graph for the appropriate altitude and range angle. 
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Table I demonstrates that this linearization 
introduces no more than 1# inaccuracy in the final 
computed values for a practical case. Special 
riot at ion seems unnecessary since the procedure 
is so simple. Any combination is readily written 
if it is remembered that each initial in-plane 
error causes four components of final error, two 
in position and two in velocity. The lateral 
errors, which are independent of the in-plane 
error, each cause two final errors, one in TRACK 
position and one in TRACK velocity. 

Lateral Position Error 

The lateral position error due to a TRACK 
error at booster burnout rates special mention 
because it illustrates an effect due entirely 
to the definition of position error used, that 
is, the type of platform assumed. The easiest 
place to view this effect is at a range angle 
of 90 degrees. It seems at first that for a 
lateral position error at booster burnout, the 
line of nodes between the plane of the actual 
ascent ellipse and the reference plane would 
lie 90 degrees away, or at apogee for this 
case. The integrated results shown in Pig. 14 
show this is not the case. 
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Pig. 16 below illustrates how this comes 
about with the central angle due to a reasonable 
position error greatly exaggerated. 




Fig. 16. Effect of Platform Definition on 
Lateral Position Error. 



Note that although no lateral velocity 
appears in the platform coordinates, with a 
platform that stays aligned with the reference 
vertical, a sidewise component of velocity, 

i 

V.^ sin U (y/R^), does show up with respect to 
the reference orbital plane. 
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CHAPTER V 
CONCLUSIONS 



If this thesis has met its objective then 
the plots showing the various error derivatives 
should clearly demonstrate the primary con- 
clusions to be drawn concerning the effect of 
target altitude and range angle on these de- 
rivatives . 

To review the results in general: 

1. For the low altitude rendezvous 
the summation of booster burnout 
velocity and injection velocity in- 
crement (XV) is only slightly great- 
er than the Hohmann value until the 
range angle is well below 90°. The 
components of follow this same 
pattern. 

2. The in-plane error derivatives 
for a 1000 n. mile target altitude 
are accurate for differentials of 

+ 60 ft/sec, + 2 milliradlans and 
12,000 feet in any direction. 



3. The lateral error derivatives are 
independent of the ln-plane error de- 
rivatives . 

4. The ln-plane partial derivatives 
are Independent of each other and there- 
fore the total error components can be 
formed by the rule for total differentials. 

5. The propulsive requirements for 
the 20,000 n. mile orbit altitude are 
shown to be so great that multi -staging 
will certainly be necessary. 

At this altitude the sensitivity 
to error is so great that special guidance 
techniques will also be required. 
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